Inconsistences of the /^-Factorization in Exclusive B-Meson Decays 



in 



O 



X 

5-H 



2 



F. Feng 1 , J.P. Ma 1 and Q. Wang 2 
1 Institute of Theoretical Physics, Academia Sinica, Beijing 100190, China 
Department of Physics and Institute of Theoretical Physics, Nanjing Normal University, Nanjing, Jiangsu 

210097, P.R.China 



Abstract 

The fcx-factorization has been widely used for exclusive decays of £>-mesons. In this factorization 
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the pertubative coefficients are extracted from scattering of off-shell partons. Because the off-shellness 
of partons the extracted pertubative coefficients in general are not gauge-invariant. We show that 
these perturbative coefficients contain gauge-dependent singularities beyond tree-level. This leads to 
that the fc^-factorization is gauge-dependent and violated in the general covariant gauge. This is not 
consistent with expectations of physics. In the factorization there exist two different wave functions for 
a £?-meson in general. We show that the two wave functions satisfy different evolution equations whose 
O solutions are used to re-sum large log's. Based on this fact, the rseummation of large log's at leading 

JL , log approximation in the widely used fc^-factorization with only one wave function is inconsistent. 

CD 

43 



Studies of exclusive £>-meson decays in experiment and theory play an important role in testing the 
standard model. For those decays involving large momentum transfers two theoretical approaches exist. 
One of them is based on the collinear factorization[lJ, in which the transverse momenta of partons in 
a hadron are integrated out and their effect at leading twist is neglected. The collinear factorization 
has been proposed for other exclusive processes for long time[2j [3]. Another one is based on the fey- 
factorization [1], where one takes the transverse momenta of partons into account at leading twist by 
meaning of wave functions. The advantage of the fey-factorization is that it may eliminate end-point 
singularities in collinear factorization [1] through resummation of large log's and some higher-twist effects 
are included. The fey-factorization, also called as pQCD approach, has been widely usedjU [6] and it is 
quite successful in phenomenology. However, the fey-factorization has not been examined so far beyond 
tree-level. 

In the fey-factorization for exclusive -B-decays the perturbative coefficients are extracted from scat- 
tering amplitudes of off-shell partons. Therefore one can not expect in general that these perturbative 
coefficients, hence the factorization, are gauge-invariant. Since the fey-factorization has not been ex- 
amined beyond tree-level for any exclusive fe?-decay, the question if the factorization is really gauge 
invariant is not answered. Using the fey-factorization for processes involving light hadrons only, e.g., the 
transition irj* — > 7, one can show that the one-loop perturbative coefficient contains gauge-dependent 
singularities [7]. The singularities are light-cone singularities and come only from the wave function of 
7r. The singularities do not appear in Feynman gauge[?l [8]. This leads to the conclusion that the fey- 
factorization is in general gauge-dependent and violated. The purpose of the current work is to examine 
this issue of exclusive decays of a l?-meson. It should be noted that in the so-called transverse momentum 
dependent factorization[9j, perturbative coefficients are extracted from scattering amplitudes of on-shell 
partons, hence they are gauge- independent. 

Another problem of the fey-factorization is related to the wave functions of 5-mesons. In the heavy 
quark limit, there exist in general two wave functions, while in the fey-factorization only one wave function 
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is employed[U [Bj. This leads to the question if the two wave functions are the same [TO]. There is no 
exact proof for the equivalence of the two wave functions. In fact, as we will show, the two wave functions 
satisfy different evolution equations. This clearly indicates that the two wave functions are not same. 
The mentioned evolution equations have the usage for the resummation of large logarithms terms in 
perturbative coefficients. Because the two wave functions and their evolution equations are different, the 
resummation of large logarithms terms in the ^-factorization with only one wave function can not be 
consistent. 

We take the semi-leptonic decay B — > tt£u£ as an example to study the above problems. We will 
use Heavy Quark Effective Theory(HQET) for the heavy 6-quark. Its field in HQET is h. The field h 
depends on the velocity v with which the S-meson moves. In the heavy quark limit, the .B-meson mass 
rriB is the same as the 6-quark mass nib, hence we have the l?-meson momentum P = mpv ~ m^v with 
v 2 = 1. The decay amplitude is determined by the matrix element, which can be decomposed into two 
form factors: 



(n(K)\q^h\B(v)) 



it 1 r 



F+ ( q 2 )(P + Kr+(m 2 B -mi)^ /•'„(, r 
F P (q 2 )P»+F K (q 2 )K». 



(1) 



In the above q is the momentum transfer q = P — K . The two form factors -Fo,+ are standard. We 
introduce two form factors Fp jc which are linear combinations of the standard two. The implication of 
Fpk will become clear later. 
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Figure 1: The leading order contribution. The thick line is for 6-quark in HQET. 

The leading order result of the fc^-factorization can be derived from the contributions in Fig.l. We 
take a frame in which the .B-meson moves in the z-direction. We will use the light-cone coordinate 
system, in which a vector o M is expressed as a M = (a + , a~, a±) = ((a + a 3 )/\/2, (a — a 3 )/\/2, a 1 , a 2 ) and 
a\ = (a 1 ) 2 + (a 2 ) 2 . We also introduce two light-cone vectors: = (0, 1, 0, 0) and = (1, 0, 0, 0). In the 
light-cone coordinate system v is given by = (v + , v _ , 0, 0). In the first step, the matrix element based 
on Fig.l can be written as: 



(n(K)\q^h\B(v)) 



1 



AN 2 
d 4 x 



(2tt)4 
—i 



dk x d 2 kx±d 4 piH^ p (p h fa) (757 + ) P7 <My, k i±) 
^■ x (0\ qi3 (x)(h v ) a (0)\B(v)), 



("Pi 



fa) 2 



r 



-v ■ fa + ie 



-ig s v 



+ 



• v \ l l ■ (K-Pl) a 

(K -pi) z + IE 



■ya , 



(2) 
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where we have already evaluated out the trivial color factors, a, /3, 7 and p are Dirac indices. We have also 
made the leading twist expansion for the quark density matrix element of ir and taken the leading order 
result with the wave function c/> n (y, ki±). Hence the momentum k% is given by ki = (0, yK~, ki±). In 
the kx factorization one sets pj = in H^g (pi, ki), the integral over p^ and z + can be then performed 
directly. One has 



(n(K)\<tfh\B(y)) 



N 2 



1 



AN 2 
dx~d 2 x± 
(2tt) 3 



dk x d 2 k 1± dp^d 2 p 1± H^ p (p 1 , h)\ p - =0 (757 + ) P7 ^(y, h±) 

e l ^--^^0\q p (x)h a (0)\B(v))\ x+=0 . (3) 



The above quark density matrix of the S-meson can be decomposed with the 16 T-matrices combined 
with the vector v^, p^j_ and or P. The appearance of the vector n or I here is because a light-cone 
direction in the density matrix is chosen. It should be noted that the two quark fields are not separated 
along a light-cone in the space-time. The decomposition can be written as: 



T a p(k + ,k±) 



dz~d 2 



(2vr) 3 
7 • v + 1 



z ± e ik+z--ik 



(0\qdz)h a (0)\B(v))\ z+=0 



7 (f> + (k + ,k±) +7 + ^ (p-(k + ,k ± ) + 



75 



(4) 



J a/3 



where the • • • stand for the terms involving 7^. The contraction of j± with H will produce perturbative 
coefficients related to these terms and the perturbative coefficients are suppressed by an extra power of 
Pi±/ m b or ku_/m b in comparison with that related to <p + or <p_. In the leading power we can neglect 
these terms. Prom the above one has two wave functions for S-meson. The exact definitions of the wave 
functions (f>+, 0_ and (f) n will be given later. With these wave functions we obtain then the leading order 
result in the /^-factorization for the form factors: 



F P (q 2 )P + = m b v + J dxdyd 2 k±d 2 p ± cf) + (x,p_ L )(j) 7r (y,k_ L )Hp(x,p ± ,y,k ± ), 
F K (q 2 )K~ = m b v~ dxdyd 2 k±d 2 p±4>-(x,p_ L )(t> 7r (y,k_ L )HK(x,p_L,y,k±), 



(5) 



with the perturbative coefficients 
H P (x,p±,y,k±) 

H K (x,p ± ,y,k±) 



N 2 - 1 
N 2 -l 



xmr 



xyrh 2 + (p± — k±) 2 ) (xrh 2 + p 2 ± ) 



+ 0(a 2 ), 



x 



p{_ 

P+ 



2N 2 
m b v~ 



y [xyrhl + (p± - k±) 2 



+ 0(a 



s)i 



~ 2 2 2 

V = m b =m b -q . 



(6) 



The integration range is given by < x < 00 and < y < 1. 

From the derivation of the leading order results several observations can be made: If we set (/>+ = </>—, 
then the above results reduces to that of [1] for the /c^-factorization of the decay in the heavy quark 
limit, where only one wave function for the .B-meson is employed. We will show later that the two wave 
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functions can not be same. In the above /cr-factorization, the partons entering the hard scattering are off- 
shell. The initial light antiquark has the momentum p^ = (xP + , 0,p±) and the fe-quark in HQET carries 
the momentum —p. The outgoing antiquark carries the momentum k^ = (0, yK~ , k±). This implies that 
one should also take these off-shell partons to calculate the higher order corrections of Hp^, as illustrated 
in [8] . A part of results of the collinear factorization can be obtained from the above results by replacing 
the wave functions with the corresponding light-cone wave functions and discard the fey-dependence in 
the perturbative coefficients and the fcy-integrals. Then we see that the so called end-point singularity 
appears in Fx with y — > 0, but not in Fp. 

The above is derived at the tree-level in Feynman gauge. In fact, the perturbative coefficients are 
already gauge-dependent at tree-level [UJ because the scattering is of off-shell partons. We will examine 
the factorization beyond the tree level. Following the above derivation we need to calculate the form 
factors and the wave functions by using the parton states with the same off-shell momenta and the 
corresponding projections for external legs indicated in Eq.(2,4) beyond the tree level. Then we can 
extract the perturbative coefficients via: 

F« = ® ® # + 4°) ® $> ® # + HP ® # , 



In this paper we will denote the n-loop contribution to a quantity A as A^ n \ From the above, one can see 

r (i) 

l P,K 



that the one-loop contribution Hp L receives contribution from the form factor and also from the wave 



functions. 

We first study the one-loop contributions from wave functions. The two wave functions of S-meson 
can be consistently defined by supplying gauge links. We introduce: 

L u (oo, z) = Pexp ^— ig s J d\u ■ G(\u + z)^J . (8) 
The wave functions can be defined by taking the limit u~ 3> n + |12l 113] : 

4>+,-(k + ,k±,0 = J ^^e lk+z --^ rk H0\q(z)Li(oo,z) l5 L u (^,0)h(0)\B(v)), 



2 = 2u-{k+f ^(u-fc) 2 
u + u 2 

with z^ = (0, z~,z±) and k + = xP + . The wave functions also depend on fj,, the renormalization scale, 
and on £ 2 . The dependence of Q 2 will be used for resummation. The limit u" S> u + should be understood 
as one neglects any contribution which is proportional to any positive power of u + . To calculate the wave 
functions consistently with the factorization, we replace the B meson with a quark pair q{p)b{—p). The 
quarks are off-shell, and p is given by p^ = (p + ,0,p±). At tree-level one simply has: 

^ ) (k + ,k 1 _)=S(k+ - P + )5\k ± -p ± ), <p^(k + ,k ± ) = 5(k + - P + )5\k ± -p ± ). (10) 

To examine the gauge dependence we will use the general covariant gauge to calculate the one-loop 
corrections to the wave functions. In this gauge the gluon propagator takes the form 

q z + ie V q z + is / 



4 




Figure 2: The one- loop contributions to wave functions. 



The Feynman gauge is obtained by taking a = 0. We will call the a-dependent corrections as gauge parts 
and denote these corrections with the sub-index a. At one- loop the corrections come from diagrams given 
in Fig.2 and Fig.3. 

As mentioned at the beginning, the one-loop corrections will contain the so called light-cone singu- 
larities. We give a detailed calculation of Fig. 2c as an example to illustrate the appearance of light-cone 
singularities. The contribution from Fig.2c to <p + and 0_ is: 



4>±(k + , kj_)\ 2c 



S(k+ -p + + q+)6 2 (k ± -p ± +q ± : 



•Tr 



(2vr) 4 

1 + 7 ■ v 
4iV, 



v • (q — p) + ie —u • q 



IE 



^ l5l ± l5 {ig s u^T a ){-ig s v"T a ) 



+ ie 



<r - a- 



+ ie 



(12) 



The integration over q + and q± can be performed with 5-functions. Performing the (/"-integral by a 
contour integral and using the pole of the heavy quark propagator at q~ = (v -p — v~q + )/v + we have for 
the gauge part: 



4>±{k + , k ± )\ 2c ,a 



k+ 



"2^ 



C F 6(-q 



v ■ p 



p' 



v + (^t(? + - p + ) + ?i) 

k± =p± -q±- 



2 ' 



(13) 



If we use it to calculate the convolution in Eq.(7) to determine Hp^ K , we have to calculate the integral 
like: 



J dk + d 2 k L f(k + M)(t>±{k + M)\2c, a K j_ dq+ J 



d 2 q±- 



<? ,rl - q±s 



l2 - 



(14) 



[{v+fql + q +( q +-p+)X 

We denote the integral in the right-hand side as /. If the function is nonzero at k + = p + and k± = pj_ 
and is approaching to zero when k + —> oo and k± —>■ oo, then we find that the integral is divergent. The 
convolution Hp^ K <S)(j)^ clearly satisfies the above requirement for the function /. Inspecting the integral 
one finds that divergence is in the region of q + — > and q± — > 0. We can regularize the divergence by 
giving the gluon a mass, or by using dimensional regularization in which the transverse momentum is in 
the space with the dimension 2 — ex,. The singularity in the integral is easily isolated by subtractions: 



f° dq + I d 2 

J — OO J 

f d " + I 

J — oo J 



q ,p±- q±) 



L[(t 



\2„2 



+ 



+f(p + ,P± 



d 2 qi_ 



f(p A 



+ g+(g H 
q + ,p± 



P +)Y [(v+y q i 



-p + )\ 



q±) 0(q 2 0± -ql)e(q+ + q+)f(p+, P± ) 



) 2 Ql 



q+p+)\ 



[{v+Yql-q+p+)\ 



f dq+ [ 



d 2 q ± 



[(v+) 2 q 2 ± -q+P+)} 2 



(15) 
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where > 0. In the above each integral is finite except the last one. The final result will not depend 
on the slicing parameter q$ and q 2 .^. The last integral is divergent. Using the integral one can identify 
the divergent part of the wave function as: 

<P±(k + ,k ± )\ 2ca = ^C F 5{k + -p + )5 2 (k ± -p L ) ( —] + finite terms. (16) 
47r \e L J 

Following our above analysis we find that the divergence appears in the momentum region of the ex- 
changed gluon whose momentum takes the patten q^ ~ (S 2 , 0(1), S, 5) with 5^0. In the following we 
will only give detailed results of singularities for the wave functions. Calculating Fig. 2a and Fig. 2b we 
also find that the contributions to 0+ contain the light-cone singularity: 

4> + (k + , k±)\2a,a = - < ^- L C F 5(k + - p + )5 2 (k± - p±) ( —] + finite terms, 

4vr \e L J 

Mk + ,k ± )\ 2b<a = ^C F 6(k + -p + )5 2 (k ± -p ± ) (J^j + finite terms. (17) 

The gauge part of the contribution from Fig. 2d to (j>+- contains an I.R. singularity regularized by the 
pole of ej = 4 — d: 

<P±(k + ,k ± )\ 2d:a = ^C F 5(k + -p + )5 2 (k ± -p ± ) + finite terms. (18) 

This I.R. singularity will be canceled by I.R. singularities from Fig. 3. 

It is straightforward to find the contributions from Fig. 2a and Fig. 2b to (ft-. It is interesting to 
note that for the contribution from Fig. 2b to 0_ we find that the gauge part is finite, while the gauge- 
independent contribution or the contribution in Feynman gauge is zero in the limit u~ 3> u~ simply 
because 7+7 • u — > 0, i.e., 

(p-(k + , k±)\2b = a [ finite terms ] . (19) 

The zero contribution from Fig. 2b in Feynman gauge leads to that the evolution equation for will be 
different than that for <p + . We come back to this later. We also find that the gauge part of 4>- from 
Fig. 2a is finite. Hence, the light-cone singularity of 0+ comes from Fig. 2a, Fig. 2b and Fig. 2c, while the 
light-cone singularity of comes only from Fig. 2c. 

Now we turn the corrections from Fig. 3. The contributions from Fig. 3a and Fig. 3d do not need to be 
considered, because the perturbative coefficients will not receive contributions from them. It is interesting 
to look at the contributions to the gauge part of <p± from Fig. 3c: 

0±(k + , fcJIaca = iag 2 C F 5(k + - p + )5 2 (k ± - p ± ) f .f q ... , ■ . "7 -— -. (20) 



(27r) 4 (q 2 + ie) 2 (v ■ (q - p) + is) 

This integral can be calculated with the standard method and the result has only an U.V. divergence. If 
we first do the (/"-integral as done for Fig. 2c, then the remaining integral is proportional to that in Eq.(15) 
by setting the function / as / = 1. In this case, unlike in Eq.(15), / will not approach to zero when 
q + — > —00 or q± — ► 00. At first look, the light-cone singularity with the momentum q^ ~ (5 2 , 0(1), S, 5) 
with 5^0 appearing in Fig. 2c will also appear here. However, in the case / = 1 one finds another light- 
cone singularity with the momentum q^ ~ (0(1), A 2 , A, A) with A — > 00. The two light-cone singularities 
cancel and it leaves only the U.V. singularity. One can also use arguments from covariance to show 
nonexistence of light-cone singularities in the above integral and in any covariant integral as shown in 
[7J. We note here that the integral in Eq.(15) with / given as Hp^ K (8) (j$ is not a covariant integral of 
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Figure 3: The one- loop corrections for the wave functions. 



q, because the range of q + with — oo < q + < is fixed at the beginning and q is fixed to be zero in 
After the U.V. subtraction the gauge parts from Fig. 3c are finite. The contributions from 
Fig.3b and Fig.3e to the gauge parts have I.R. singularities. The results from Fig.3c, Fig.3b and Fig.3e 
are: 

4>±{k + , A;_l)| 3CiQ = a [ finite terms ] , 

<t>±{k + M)hb,a = 0±(k + ,k ± )\ 3e , a = -^C F 5(k + -p + )5\k ± -p ± )(^j+ finite terms. (21) 

Similarly, there is no contribution from Fig.3f to (/)- in Feynman gauge. The gauge part of (j)+- has 
only U.V. divergence. After the U.V. subtraction we have: 

4>+{k + , k±)\ 3 f, a = a[ finite terms ] , 
4>-(k + ,k ± )\ 3f = a[ finite terms ] . (22) 

Adding everything together we find that the wave functions contain the light-cone singularity: 



i 1} 



(k + , k±) = ^C F 5(k + - p + )5 2 (k ± - p ± ) (— ) + finite terms. 
4vr \eL, 



(23) 



The I.R. singularity from Fig. 2d is canceled by that from Fig. 3b and Fig.3e. 

We turn to the wave function of ir. For ir with the momentum = (0, K~ , 0, 0) one can define the 
wave function with the gauge link along the direction u with = (u + , u~ , 0, 0). The wave function is 
defined in the limit u + S> u~ as 

f^^^^ + -^^(0|g(0)4(oo.O) 7 - 75 L,(oo.z) (? (^)| 7 r(AO)| z - = o, 

(24) 



k~ _ xK~ , ..«*■*? 



ti- 



ll 2 



We take a quark pair to replace the pion and calculate the wave function. The quark q of the pair carries 
the momentum k\ and the antiquark q carries the momentum K — k±. k\ is given as k^ = (0, fcuj 
with ki = xqK~. At leading order we have: 



<t>W(x,k ± ) = 6(k~ - kY)5(k ± - k 1± ). 



(25) 
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At one-loop level the wave function receives corrections from the same diagrams in Fig. 2 and Fig. 3. by 
replacing the 6-quark line with the light quark line. Calculating them in the general covariant gauge, 
we also find the light-cone singularities. From our previous results in [7] converted with the dimensional 
regularization we have the gauge parts containing the light-cone singularity: 



\{k~, k±)\ 2a ,a 
^(Ar ,k±)\ 2c>a 



Cpb{k — k l )5 (k± — ki±) ( — ) + finite terms, 
47r ~ \e L , 



n n , ,., _ 2-"' ■ ' / 2 



-| — - — Cp5(k — k 1 )S (k± — kx±) — + finite terms, 

H — —^-CpS(k~ — ki)5 2 {k± — k\x) [ — ) + finite terms, 

aa, „ _ , _x r o/r j* \ ( 2 , , . 

1 + hmte terms. 



-C F 5(k- -kTW(k ± -k 1± ) 
4vr V< /. 



(26) 



To extract the one- loop perturbative coefficients one needs to use the same quark pairs, used to 
calculate the wave functions, to calculate the form factors. At one-loop level, the form factors will 
not have such light-cone singularities, as we meet in wave functions. The reason for this is similar to 
the case of Fig. 3c, where the loop integral is covariant. With the argument from the covariance one 
can show that scattering amplitudes, i.e., the form factors, do not have the light-cone singularity in 
the general covariant gauge. Details can be found in [7j. This is in agreement with the expectation 
that any scattering amplitude has only U.V. divergences, I.R. divergences and collinear divergences. 
Actually, scattering amplitudes with off-shell partons will not have any soft divergence like collinear- 
and I.R. singularity, because they are all regularized by the off-shellness of partons. Hence the one loop 
perturbative coefficients will have the light-cone singularities: 



Hp\x,p_L,y,k±) 
Hx\x,p±,y,k±) 



~-H { p\x,p±,y, k±) 
■^-HP(x,p ± ,y,k±) 



1 )-> 



+ finite terms, 
+ finite terms. 



(27) 



It should be noted that the singular contributions depend on gauges. Therefore, the /c^-factorization is 
gauge dependent. Also because it is singular, the /cp-factorization can not hold beyond tree-level in the 
general covariant gauge. 

Now we turn to the ^-dependence of <f>±. This dependence can be used to do resummation of large 
log's. The dependence of <f>+ can be found in [13], where it is determined with on-shell quark pair and 
in Feynman gauge. From our calculation one can verify that the dependence is gauge independent. The 
result for d> + is 13 : 



4as 
' 3vr 



In 



C 2 6 2 e 2 7 -1 



2a s n 2 e 
in — 7T- 

3vr C 2 



2a, 



H 2 e 



K(jji,b) + GQi,Q--^-)n^ 



(f>+(k + ,b, /i,0 
(f) + (k + ,b,fi,(), 



(28) 



it should be noted that the kernel is not exactly the famous factor K + G, the last factor comes because 
we used HQET for the heavy quark. In the above the wave function is defined in the impact space as: 



0±(fc + ,6,/x,C) 



d 2 k ± e ik ^ 



b <t>±{k + , fc_L,/x,C). 



(29) 
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The (-dependence is introduced by those diagrams in Fig. 2 and Fig. 3, where the exchanged gluon is 
emitted or absorbed by the gauge link. Only from Fig. 2b and Fig.3f, <^>+ receives contributions which 
are proportional to ln 2 ( 2 . This results in that Eq.(28) contains a single log of ( 2 . This dependence is 
useful for resummation of Sudkov double log's. As we already noticed before that Fig. 2b and Fig.3f do 
not contribute to 0_ in Feynman gauge. Except these two diagrams, other diagrams contribute to the 
same (-dependence to 0_ as to <p + . Using our previous results in [13] we find the evolution equation for 

0_(fc+,6,/2,C). (30) 

From the above, the (-dependence of (j>- is different than that of Therefore we in general can not 
have 0+ = </>-• The (-dependence of (/>_ is a constant of (, it can only be used to re-sum single log's. 
By using the solutions of the above evolution equations one can resum possible log terms in perturbative 
coefficients into an exponential factor S. Since the evolution equations are different, the exponential 
factor S related to 0+ or the form factor Fp is different than that related to (ft- or the form factor Fx- 
This implies that the resummation in the A>r-factorization with only one wave function in [U E] can 
not be consistent. 

To our conclusion: The fcr-factorization has been widely used for exclusive -B-decays and it is quite 
successful in phenomenology. However, its successfulness does not imply theoretical consistence of the 
factorization. In the /^-factorization of exclusive -B-decays the perturbative coefficients are extracted 
from scattering amplitudes of off-shell partons. Because the scattering amplitudes are gauge-dependent, 
one can not expect in general that the extracted perturbative coefficients are gauge invariant. We have 
taken the semi-leptonic decay B — * ir£i>£ as an example to show that the perturbative coefficients are 
indeed gauge-dependent. They contain light-cone singularities in the general covariant gauge. The 
singularities disappear in Feynman gauge. The singular contributions come only from wave functions. 
This leads to the conclusion that the fc^-factorization of exclusive -B-decays is gauge dependent and 
violated in the general covariant gauge because of the singularities. A gauge-dependent factorization 
is generally unacceptable since it deliveries physical predictions which are gauge-dependent. The kx- 
factorization employs only one wave function to describe the nonperturbative property of a -B-meson. 
In fact there are two wave functions. The two wave functions are different and they satisfy different 
evolution equations which can be used to re-sum large log's. Because the difference of the two wave 
functions and their evolution equations, the resummation of large log's at leading log approximation in 
the ^-factorization with only one wave function is inconsistent. 
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